Nanoscale surface relaxation of a membrane stack 
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Recent measurements of the short-wavelength (~ 1-100 nm) fluctuations in stacks of lipid mem- 
branes have revealed two distinct relaxations: a fast one (decay rate of ~ 0.1 ns~^), which fits the 
known baroclinic mode of bulk lamellar phases, and a slower one (~ 1-10 fJ-s~^) of unknown origin. 
We show that the latter is accounted for by an overdamped capillary mode, depending on the surface 
tension of the stack and its anisotropic viscosity. We thereby demonstrate how the dynamic surface 
tension of membrane stacks could be extracted from such measurements. 
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Self-assembled stacks of membranes are encountered in 
various industrial and biological systems. They consist 
of parallel bilayers of amphiphilic molecules separated by 
microscopic layers of solvent — a structure with the sym- 
metry of a smectic A liquid crystal [3] . Such stacks form 
lyotropic lamellar phases 0] , on which many cleaning and 
cosmetic products are based. Lamellar bodies arc found 
also in the lung 0] and as multilayer vesicles ( "onions" ) 
Membrane stacks made of phospholipids have been 
widely used to study properties of biological membranes, 
whereby the large number of identical, equally spaced 
membranes helps enhance the signal and allows the study 
of membrane-membrane interactions (e.g., Q). 

The elasticity of membrane stacks is equivalent to that 
of single-component (thermotropic) smectics [l| and has 
been extensively studied. The elastic moduli of the stack 
can be extracted from its equilibrium fluctuations using, 
e.g., x-ray line shape analysis [6j|. By contrast, the hydro- 
dynamics of membrane stacks 0, Q , because of their two 
micro-phase-separated components, differs from that of 
thermotropic smectics An additional hydrodynamic 
mode appears — the baroclinic (slip) mode — along with 
a unique dissipation mechanism, in which the membranes 
and solvent layers develop different average velocities Q ■ 
Experimental studies of hydrodynamic modes in mem- 
brane stacks have been rather scarce, the prevalent tech- 
nique being dynamic light scattering 8], whose spatial 
resolution is limited by the wavelength of light. 

In a recent experiment using neutron spin-echo spec- 
trometry, Rheinstadter, Haufiler, and Salditt (RHS) have 
provided a first look at the relaxation of membrane stacks 
at short wavelength (1-100 nm) and short time (1-10"^ 
ns) [l3|- Their system consisted of several thousands 
of dimyristoylphosphatidylcholine (DMPC) phospholipid 
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bilayers, self-assembled into a stack of d ^ 5 nm period- 
icity. The system was studied at temperatures above 
and below the lipid melting point, corresponding to fluid 
and gel-like membranes, respectively. In both cases the 
measured dynamics consisted of two distinct exponential 
relaxations. The dispersion relation of the faster relax- 
ation (decay rate of 0.1 ns""'^) could be well fitted in 
the fluid- membrane case to that of the baroclinic mode of 
a bulk lamellar phase [lH, while the slower mode (decay 
rate of ~ 1-10 /is~^) was left unexplained. We demon- 
strate below that this slower relaxation is well accounted 
for by a surface mode, i.e., a perturbation which is local- 
ized within a finite penetration depth from the surface of 
the stack. 

In a recent publication [l2| we have addressed the sur- 
face dynamics of membrane stacks, highlighting the qual- 
itative differences from the surface dynamics of both sim- 
ple liquids and thermotropic smectics [13]. These differ- 
ences arise from the slip dissipation mechanism, which is 
absent in simple liquids and thermotropic smectics but is 
usually dominant in lyotropic lamellar phases. Although 
the formulation in Ref. is general, its analysis is fo- 
cused on a very different domain (larger wavelengths and 
slower rates) from that sampled by RHS. In that do- 
main the slip dissipation dominates and, consequently, 
the surface relaxation is governed by an overdamped dif- 
fusive mode, whose decay rate F increases quadratically 
with the wavevector q. In this Brief Report we present 
a slight adaptation of that theory for a large-q, high-F 
regime such as that of RHS. 

The general surface dynamics of membrane stacks is 
quite complex, depending on several restoring and dis- 
sipation mechanisms [l^. Three moduli are associated 
with the restoring forces: the compression modulus B, 
bending modulus K , and surface tension 7. Viscous dis- 
sipation is characterized (in the limit of incompressible 
flow) by three viscosity coefficients 7], denoted rjM, iJt, 
and 77V ■ The coefficient 77M, associated with differences 
in the lateral velocity across layers (sliding viscosity), is 
much smaller than the other two, which correspond to 
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the viscous response to deformations of the hpid mem- 
branes. We use the parameter 6 = 2(77t + ?7v)/w to 
characterize this viscosity anisotropy; it is typically of 
order 10^-10'^ 0, [13] ■ The aforementioned slip motion 
requires another transport coefficient A* — d'^/{l2r]o), 
where ijo is the viscosity of the solvent (water) layer. 

In view of this richness it is helpful to begin by iden- 
tifying the dominant contributions to the slower mode 
of Ref. First, for the typical parameters of that 

case — q ^ 10~^ nm^^, F ~ 1 fJ,s~^, ijm ~ 10~^ Pa s, 
and mass density p ^ 1 g/cm'^ — one gets a negligible 
Reynolds number, Re ~ pT/(riyiq^) ~ 10~^, implying 
that inertial modes are irrelevant in the current case. 
Second, to determine the dominant dissipation mecha- 
nism one should compare the friction due to slip, 
{v being a characteristic relative velocity), with that due 
to viscous stresses, rnyq^v, i.e., the dimensionless pa- 
rameter S — {rjMlJ-q^)~^ is to be compared with Q 
We find 5 ~ 10 ^ 0. Thus, unlike the mode focused 
on in Ref. in the current large-g case viscous dissi- 
pation is dominant. Finally, the relative importance of 
the three restoring mechanisms depends not only on the 
surface perturbation wavevector q but also on its pen- 
etration depth . Since the value of a is unknown a 
priori, all three mechanisms should be considered in prin- 
ciple. However, to keep the analysis as simple as possible 
we shall assume that the surface tension is the dominant 
factor. This ansatz is motivated by the experimental fact 
that the rate of RHS's slower mode is linear in \q\ at small 
q (see Fig. [Ij ; the way to get such a linear overdamped 
dispersion relation is to balance a surface tension stress 
against a viscous one, ^q^u ~ rjqTu (u being the ampli- 
tude of the surface deformation). We will return to the 
consistency of this assumption later on. 

The continuum theory formulated in Ref. ^12i| is valid 
for wavelengths much larger than the inter-membrane 
spacing, qd <g; 1. RHS's experiment, however, samples 
the range 0.1 < qd < 4. To obtain an extrapolation 
of the analysis to large q we introduce one last modifi- 
cation to the theory — the distance z from the surface 
into the stack is discretized, z — > —dn {n = 0, 1, 2, . . .), 
turning the differential equations of Ref. [l2| into finite- 
difference ones (similar to the analysis of high-g acoustic 
modes in a crystal). The lateral position x parallel to 
the membranes is kept continuous, and we consider, for 
simplicity, a surface perturbation which is uniform in the 
second lateral direction y. 

Within these assumptions Eq. (11) of Ref. yields 
the following surface mode for the vertical displacements 
of the membranes, u„(a;,t): 

Un = (C+e-"+''"-KC_e-"-''")e*'^-^* 
2 



(1) 



For sufficiently small q (qd <^ 0^^/^) the spatial decay 
coefficients are a± ~ i.e., the mode contains 

two terms of disparate penetration depths, aZ^ 3> aT^. 



(A qualitatively similar result was obtained for the sur- 
face mode analyzed in Ref. [l3], yet in the current case 
the origin of the two differing penetration depths is the 
large viscosity anisotropy rather than the strong slip dis- 
sipation.) In the other limit of qd ^ Q^^^, as expected, 
both contributions become localized within a distance of 
order d from the surface, a± ~ (2/(i) ln(0^^/'^|(7|c?). 

The dispersion relation T{q) is set by the boundary 
conditions for the stress tensor at the stack surface, as 
summarized in Eq. (13) of Ref. [l^. Substituting in that 
equation the expressions for a± obtained above, we get, 
within the same approximations, 
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Qriud 
sinh^i (-Q^/^\q\d 



sinh^i Qe-i/2|^|rf^ 



•(2) 



Equation ^ is the main result of our current analysis. 
For large wavelengths this dispersion relation becomes 



2[ryM(r/T+w)/2]i/2 



(3) 



Equation ([3]) is equivalent to the dispersion relation of 
an overdamped capillary mode at the surface of a simple 

liquid having effective viscosity r^eff = [f?M('7T+w)/2]^^^- 
In the opposite, short-wavelength limit we get 



T{qd > e^/') 



H\q\d). 



(4) 



In this quasi-two-dimensional limit the dependence on 
the smaller (sliding) viscosity, 77M, disappears, and an 
effective two-dimensional viscosity emerges, 7^20 = (vt + 

Figure [T] shows fits of the dispersion relations for the 
slower mode, as measured by RHS, to Eq. ^ pjj]. (The 
measurements for q > 0.5 nm~^ are considered less reli- 
able due to scattering by defects in the stack [13 •) The 
stack periodicity was measured as d = 5.4 and 5.6 nm 
at temperatures T — 30°C (fiuid membranes) and 19°C 
(gel- like membranes), respectively [l0[- The value of the 
shding viscosity at 30°C, riu = 0.016 Pa s, was indepen- 
dently found from a fit of the faster mode [l3| . We are 
thus left with two fitting parameters in Eq. O and 7. 
For the fluid-membrane case we find O = 110 and 7 — 5.4 
mN/m. It should be stressed that having two free param- 
eters does not allow for accurate determination of both, 
and these values should be regarded merely as rough es- 
timates. Nonetheless, the fitted values are of the correct 
scale. The value for 8 implies r?T,v ~ 0f?M 1 Pa s, 
i.e., a viscosity 3 orders of magnitude larger than that 
of water, which matches the typical effective viscosity of 
lipid membranes [T7| . It also implies an effective two- 
dimensional viscosity 7720 ~ VT.yd ^ 10~^-10~^ Pa s m, 
which agrees well with measurements of the surface vis- 
cosity of fiuid DMPC membranes p^. 

The applicability of the theory to stacks of solid, gel- 
like membranes should be questioned, as such stacks have 
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FIG. 1: Dispersion relations for the slower relaxation mode of 
stacks of DMPC lipid membranes at 30° C (circles) and 19° C 
(squares). (Data taken from Ref. P^].) The solid lines are 
fits to Eq. © with d = 5.4 nm, r?M = 0.016 Pa s, 9 = 110, 
and 7 = 5.4 mN/m (lower curve); d = 5.6 nm, tjm = 0.016 
Pa s, B = 350, and 7 = 28 mN/m (upper curve). The values 
of d and t^m are taken from Ref. lO]; O and 7 are fitting 
parameters. 

additional intra-membrane elasticity. The same concern, 
in fact, should be raised regarding the fluid-membrane 
case as well, since at the high frequencies considered 
here the individual membranes are expected to have a 
viscoelastic response. The fits obtained in Fig. [T] (in par- 
ticular, the linear behavior for small q) suggest, however, 
that these additional restoring forces are negligible com- 
pared to the surface tension and do not affect the sur- 
face relaxation. The fit for T = 19° C yields significantly 
larger values for both the viscosity anisotropy and the 
surface tension, 8 — 350 and 7 = 28 mN/m, which is 
the expected trend for stiffer membranes [19[ . (In the fit 
we have assumed that the sliding viscosity rjyi does not 
change much with temperature.) 

The elasticity of membrane stacks gives rise to an ef- 
fective static surface tension, 701 — {KBY^^ [H. [20|. [21| . 
The values of K and B in the fluid-membrane state were 
extracted by RHS from the fit of the faster mode as 
K ~ 1.15 X 10"^^ N (corresponding to a membrane bend- 
ing modulus K ^ Kd ~ 14.8 /cbT) and B ~ 1.08 x 10^ 
Pa. This yields 7ci — 11.1 mN/m, which is of the same 



order of magnitude as the high-frequency surface tension 
fitted above [2^. 

We now return to check the self-consistency of our as- 
sumptions. First, for a mode to be a surface one its 
penetration depth must be smaller than the total thick- 
ness of the sample. The penetration depth found in Eq. 
dH) is aZ^ < e^/2/g < 10^ nm, which is at least 1-2 
orders of magnitude smaller than the thickness of RHS's 
films (~ 10 fim). Second, for surface tension to be the 
dominant restoring force, one should have 7 > 7ei. This 
condition can be obtained rigorously (23| but is also real- 
ized upon demanding that the stress arising from surface 
tension, jq^au, be larger than both the compression one, 
Ba^u, and the bending one, Kq^u. As described above, 
we actually have 7 ~ 701 and, thus, the assumption can 
be only marginally fulfilled. Moreover, the omission of 
the bending terms requires also that Kq^ j^yy be smaller 
than r jj^, which is satisfied only for the lowest end of 
the sampled q range, q ^ 0.1 nm^^. The apparent suc- 
cess of the simplified theory over the extended q range 
(Fig. [T|), therefore, is somewhat surprising. We note that 
the stacks of RHS are densely packed. The thickness of 
a DMPC bilayer at 30°C is 4.5 nm [2|], implying that 
the solvent layers in-between membranes are only 1 nm 
thick. For such density and high-g surface perturbations 
the stack might not follow the usual description of linear 
smectic elasticity but respond merely as an anisotropic 
viscous liquid with surface tension. 

In summary, the relaxation of nanoscale fluctuations 
in finite membrane stacks seems to occur via two dis- 
tinct overdamped modes — a bulk baroclinic mode and 
a slower surface mode. The dispersion relation of the sur- 
face mode provides access to the dynamic surface tension 
of the stack, which should be hard to measure otherwise. 
Supplementing such an experiment with measurements 
at larger wavelengths (e.g., using dynamic light scatter- 
ing), yielding a value for 0, may allow the accurate ex- 
traction of the dynamic surface tension. 
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